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on behalf of the ARGO-YBJ coUaboration 

Important informations on the origin and the propagation mechanisms of cosmic rays may be 
provided by the measurement of the anisotropies of their arrival direction. In this paper the obser- 
vation of anisotropy regions at different angular scales is reported. In particular, the observation of 
a possible anisotropy on scales between ~10° and ~30° may be a key-detection for speculations on 
the presence of unknown features of the magnetic fields the charged cosmic rays propagate through, 
as well as to potential contributions of nearby sources to the total flux of cosmic rays. Evidence of 
new weaker few-degree excesses throughout the sky region 195° < R.A. < 315° is also reported. 



Introduction 

As the most part of cosmic rays (CRs) are charged 
nuclei, their arrival direction is deflected and made 
isotropic by the action of galactic magnetic field 
(GMF) that they propagate through before reaching 
the Earth atmosphere. In such a field, the gyroradius 
of CRs is given by ra.u. — 100i?Tv, where ra.u. is in 
astronomic units and Rty is in TeraVolt. It must be 
taken as a reference value, because the GMF is the 
superposition of regular field lines and chaotic contri- 
butions, the strength of both them being still under 
debate. Data available today gives for the local total 
intensity the value 5 = 2-^4 /iG. 

Actually, deviations from the isotropy are expected 
to occurr as a consequence of the particular realiza- 
tion of the random distribution of cosmic ray sources 
and magnetic field lines in the galaxy [l|. In this 
framework, the amplitude of the anisotropy is propor- 
tional to the rigidity, that is why it is expected to be 
seen at high energies (above few hundreds TeV). How- 
ever, different experiments 0-0| observed an energy- 
dependent "large scale" anisotropy in the sidereal 
time frame, well below that threshold. The ampli- 
tude is about 10~* - 10"'^, suggesting the existence 
of two distint broad regions, one showing an excess of 
CRs (called "tail-in"), distributed around 40° to 90° 
in Right Ascension (R.A.). The other a deficit (the 
"loss cone"), distributed around 150° to 240° in R.A.. 
The origin of these anisotropies is still unknown. Some 
authors claim that the observations may be due to a 
combined effect of the regular and turbolent GMF Q , 
or to local uni- and bi-dimensional inflows [o!]. Other 
studies suggest that it can be explained within the 
diffusion approximation taking into account the role 
of the few most nearby and recent sources [l|, [l3| • 

Easy to understand, more beamed the anisotropics 
and lower their energy, more difficult to fit the stan- 
dard model of CRs and GMF to experimental results. 
That is why the evidence of the existence of a medium 
angular scale anisotropy contained in the tail-in re- 
gion by the Tibet AS7 [Til and Milagro [12} collabo- 



rations in rcent years was rather surprising. Similar 
small scale anisotropies has been recently claimed to 
be observed by the Icecube experiment in the South- 
ern hemisphere Q- So far, no theory of CRs in the 
Galaxy exists which is able to explain few degrees 
anisotropies in the rigidity region 1-10 TV leaving the 
standard model of CRs and that of the local GMF 
unchanged at the same time. 

From the experimental viewpoint, observing 
anisotropy effects at the level of lO""* with an air 
shower array is a difficult job, because of the intrin- 
sic difficulties that this kind of apparatus has to cope 
with in estimating the exposure. 

Finally, the observation of a possible small angular 
scale anisotropy region contained inside a larger one 
rely on the capability for suppressing the anisotropic 
structures at larger scales without, at the same time, 
introducing effects of the analysis on smaller scales. 

In this paper the observation of CR anisotropy at 
different angular scales with ARGO-YBJ is reported 
as a function of the primary energy. 



I. THE ARGO-YBJ EXPERIMENT 

The ARGO-YBJ experiment, located at the Yang- 
BaJing Cosmic Ray Laboratory (Tibet, P.R. China, 
4300 m a.s.l., 606 g/cm^), is an air shower array able 
to detect the cosmic radiation at an energy threshold 
of a few hundred GeV. The full detector is in stable 
data taking since November 2007 with a duty cycle 
greater than 85%. The trigger rate at the threshold 
is 3.6 kHz. The detector characteristics and perfor- 
mance are described in [13l] . 



II. DATA ANALYSIS AND RESULTS 

In order to study the anisotropy at different angu- 
lar scales the isotropic background of CRs has been 
estimated with two methods: the equi-zenith angle 
method [l3] and the direct integration method [Tsj . 
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The equi-zenith angle method, used to study the 
large scale anisotropy, is able to eliminate various spu- 
rious effects caused by instrumental and environmen- 
tal variations, such as changes in pressure and temper- 
ature that are hard to control and tend to introduce 
systematic errors in the measurement. The method 
uses data coming from all the angular scales, so that 
potential small structures are not separated from the 
underlying large scale modulation. 

The direct integration method, based on time- 
average, rely on the assumption that the local distri- 
bution of the incoming CRs is slowly varying and the 
time-averaged signal may be used as a good estimation 
of the background content. Time-averaging methods 
act effectively as a high-pass filter, not allowing to 
inspect features larger than the time over which the 
background is computed (i.e., 15°/hourxAi in R.A.). 
The time interval used to compute the average spans 
Ai= 3 hours and makes us confident the results are 
reliable for structures up to «35° wide. 



A. Large Scale Anisotropy 

The observation of the CR large scale anisotropy by 
ARGO-YBJ is shown in the figure [1] at different pri- 
mary energy up to about 25 TeV. The data used in this 
analysis was collected by ARGO-YBJ from 2008 Jan- 
uary to 2009 December with a reconstructed zenith 
angle < 45°. The so-called 'tail-in' and 'loss-cone' re- 
gions, correlated to an enhancement and a deficit of 
CRs, are clearly visible with a statistical significance 
greater than 20 s.d.. The tail- in broad structure ap- 
pears to dissolve to smaller angular scale spots with 
increasing energy. It should be stressed that the ener- 
gies reported in the figure[T]refer to the median energy 
of all nuclei triggering the experiment. 

To quantify the scale of the anisotropy we studied 
the 1-D R.A. projections integrating the sky maps in- 
side a declination band given by the field of view of 
the detector. Therefore, we fitted the R.A. profiles 
with the first two harmonics. The resulting ampli- 
tude of the first harmonic is plotted in the right plot 
of figure [2] where is compared to other measurements 
as a function of the energy. The ARGO-YBJ results 
are in agreement with other experiments suggesting 
a decrease of the anisotropy first harmonic amplitude 
with increasing energy. 




FIG. 1: Large scale CR anisotropy observed by ARGO- 
YBJ as a function of the energy. The color scale gives the 
relative CR intensity. 
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FIG. 2: Amplitude of the first harmonic as a function of 
the energy, compared to other measurements. 



B. Intermediate Scale Anisotropy 

The figure [3] shows the ARGO-YBJ sky map in 
equatorial coordinates. The analysis refers to events 
collected from November 2007 to May 2011 after the 
following selections: (1) >25 shower particles on the 
detector; (2) zenith angle of the reconstructed showers 



<50°. The triggering showers that passed the selec- 
tion were about 2-10^^. The zenith cut selects the 
declination region 5 ~ -20°-^ 80°. According to the 
simulation, the median energy of the isotropic cosmic 
ray proton flux is «1.8 TeV (mode energy «0.7 
TeV). 
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FIG. 3: Intermediate scale CR anisotropy observed by ARGO-YBJ. The color scale gives the statistical significance of 
the observation in standard deviations. 
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FIG. 4: Size spectrum of the regions 1 and 2. The vertical 
axis represents the ratio between the events collected. The 
upper scale shows the corresponding proton median energy 
(see text). 

The most evident features are observed by ARGO- 
YBJ around the positions a ~ 120°, 5 ~ 40° and 
a ^ 60°, 5 ^ -5°, positionally coincident with the re- 
gions detected by Milagro [l^]. These regions, named 
"region 1" and "region 2", are observed with a sta- 
tistical significance of about 14 s.d.. The deficit re- 
gions parallel to the excesses are due to a known ef- 
fect of the analysis, that uses also the excess events 
to evaluate the background, artificially increasing the 
background. On the left side of the sky map, several 
possible new extended features are visible, though less 
intense than those aforementioned. 

The area 195° < R.A. < 315° seems to be full of 



few-degree excesses not compatible with random fluc- 
tuations (the statistical significance is more than 6 s.d. 
post-trial). The observation of these structures is re- 
ported here for the first time and together with that of 
regions 1 and 2 it may open the way to an interesting 
study of the TeV CR sky. 

To figure out the energy spectrum of the excesses, 
data have been divided into five independent shower 
multiplicity sets. The number of events collected 
within each region are computed for the event map 
as well as for the background one. The ratio of these 
quantities is computed for each multiplicity interval. 
The result is shown in the figure 2] Region 1 seems 
to have spectrum harder than isotropic CRs and a 
cutoff around 600 shower particles (proton median 
energy Ep° = 8 TeV). On the other hand, the excess 
hosted in region 2 is less intense and seems to have 
a spectrum more similar to that of isotropic cosmic 
rays. The steepening from 100 shower particles on 
(Ep° — 2 TeV) is likely related to efficiency effects. 
Further studies are on the way. 

The figure [5] reports the amplitude of the regionl 
and 2 anisotropics as a function of the UT. Time unit 
is the Modified Julian Date. It can be appreciated 
that no time evolution of the anisotropics are there, 
nor evidence of correlation of the emission. 



III. CONCLUSIONS 

This paper reports the observation of CR anisotropy 
at different angular scales with ARGO-YBJ, as a 
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FIG. 5: Region 1 (upper plot) and 2 (lower plot) relative 
intensity as a function of UT. 



function of the primary energy. The large scale CR 
anisotropy has been clearly observed up to about 25 
TeV. Evidence of existence of different few-degree ex- 
cesses in the Northern sky (the strongest ones posi- 
tionally coincident with the regions detected by Mila- 
gro in 2008) is reported. The time distribution of the 



phenomenon has been showed too. A discussion of the 
systematic effects which may come from imprecisions 
in estimating the reference level is outlined. 

In fact, a wrong estimation affect the significance 
and the relative intensity sky maps, even creating ar- 
tifacts (i.e. fake excesses or deficit regions). Drifts in 
detector operating conditions or atmospheric effects 
on the air shower development are quite hard to be 
modeled, then difficult to be accounted for to suffi- 
cient accuracy. As anisotropics of the order 10"'' are 
looked for, operating conditions must be known down 
to this level, all across the field of view and during all 
the acquisition time. 

Background methods applied in this analysis 
demonstrated to be effective for analysis of diffuse 
from even weaker regions (e.g. see fi^). 

Given the importance of the topic, a joint analy- 
sis of concurrent data recorded by different experi- 
ments in both hemispheres, as well as a correlation 
with other observables like the interstellar energetic 
neutral atoms distribution [13, [ill ' should be a high 
priority to clarify the observations. 
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